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ABSTRACT 

Edge Cloud 2 (EC2) is a molecular cloud, about 35 pc in size, with one of the largest galac- 
tocentric distances known to exist in the Milky Way. We present observations of a peak CO 
emission region in the cloud and use these to determine its physical characteristics. We calculate 
a gas temperature of 20 K and a density of ri(H2) ~ lO^cm"^. Based on our CO maps, we 
estimate the mass of EC2 at around 10'* Mq and continuum observations suggest a dust-to-gas 
mass ratio as low as 0.001. Chemical models have been developed to reproduce the abundances 
in EC2 and they indicate that: heavy element abundances may be reduced by a factor of five 
relative to the solar neighbourhood (similar to dwarf irregular galaxies and damped Lyman alpha 
systems); very low extinction {Ay < 4 mag) due to a very low dust-to-gas ratio; an enhanced cos- 
mic ray ionisation rate; and a higher UV field compared to local interstellar values. The reduced 
abundances may be attributed to the low level of star formation in this region and are probably 
also related to the continuing infall of primordial (or low metallicity) halo gas since the Milky 
Way formed. Finally, we note that shocks from the old supernova remnant GSH 138-01-94 may 
have determined the morphology and dynamics of EC2. 

Subject headings: ISM: clouds, dust, molecules - Galaxy: edge, far-outer - clouds: individual (Edge 
Cloud 2) - astrochemistry: models - radio lines: ISM 

1. Introduction 



inow at NRAO, Green Bank, WV. 

2now at School of Mathematics and Physics, Queen's Observations of CO emission at large galacto- 

Universlty Belfast. Centric distances have detected a number of molec- 

^now at Department of Astronomy, Yale University, ular clouds, among which is Edge Cloud 2 (here- 
New Haven, CT. j^f^gj. with an estimated kinematic galacto- 



centric distance of 28kpc (jPigel et al.lll994[ ). 
Although this kinematic distance for the far-outer 
Galaxy should be treated with caution, EC2 ap- 
peared to be up to 6 kpc further away than Edge 
Cloud 1 (ECl), t he next mos t distant molecular 
cloud detected by iDigel et aL , and much further 
than the extent of the optical disk of the Milky 
Way ~ 15-19 kpc (jFich et alJ Il989l : iRobin et aL 
19921 ). This is almost as far as the most distant 



Hi, detected .a-t ^ 30 kpc (iKulkarni et al.lll982f) . 



More recently iLevine et al. ( 2006 ) have traced the 
spiral structure in the southern half of the Galaxy 
out to at least 25 kpc, implying a minimum ra- 
dius for the gas disk. EC2 has an effective radius 
of ~20 pc and is situated ^3 50 pc below the dis - 



tant warped Galactic plane (jPigel et al.l Il996b[ ). 



The CO luminosity of EC2 is at least a factor of 2 
larger than those of the other 10 clouds detected 
at large distances, and is compa rable to that of th e 
Taurus Giant Molecular Cloud (jPigel et al.lll994) . 
EC2 is also the only e dge cloud detected in the 



high-density tracer CS (jPigel et al.lll996al) . The 



CO maps of EC2 show that it has sub-structure. 
It is in a region of extremely low gas pressure and 
very small spiral arm perturbation. 

EC2 was also found to have an associated 
Hll region excited by an early B star (MRl) 
( de Geus et al. 19931 ) that appears to have trig- 
gered star formation. Studies of metallicity as 
a function of galactocentric distance have shown 
that there is a galactic metallicit y gradient. Spec- 
tra of MRl dSmartt et all Il996 l) indicate signifi- 
cant metal depletion - with elemental abundances 
reduced on average by some 0.5 dex - and that its 
galactocentric distance is more likely to be in the 
range 15-19 kpc (placing it at the limit of t he op - 
tical disk). Subsequently, iRoUeston et al. I (|2000l) 
have calculated metal dep letions of about five for 
C, N, and O. Kobayashi fc Tokunaga (2000) used 
NIR observations to argue that MRl has triggered 
the format i on of young stellar objects in EC2. 
Snell et al.l (l2002l) argue that it is the most distant 
cloud in the Milky Way, with evidence for massive 
star formation, alth ough this distinction is c laime d 
for WB89-789 by iBrand fc WouterlootI (|2007l ). 
More recently, two embedded yo ung star clusters 



have been associated with EC2 (Kobavashi et al 



20051 ). lYasui et al.l (|2006f) identified 52 members 
in one of the clusters, sug gesting that t he clu ster 
is a T Tauri association. IStil fc IrwinI (|200lh as- 



sociate EC2 with the approaching side of the H I 
shell from supernova remnant (SNR) GSH 138- 
01-94, thereby setting the galactocentric distance 
R ^ 23.6 kpc for EC2. Because of the uncertainty 
of R for EC2, we have used two heliocentric dis- 
tances, 14 and 20 kpc {R = 22 and 28 kpc) for any 
distance sensitive calculations in this work. 

Current models of galacto-chemical evolution 
include time and spatial variations in the halo gas 
infall and star formation rate whereby the Galactic 
halo, bulge, and thick disk formed first, separately 
from the thin disk, in two infall episodes. Chem- 
ical abundances of the interstellar medium and 
their radial variation across galactic disks provide 
a fundamental set of constraints for theories of 
disk formation and evolution. The most accepted 
mechanism to explain the existence of abundance 
gradien ts in disk galaxies is the so-ca lled 'biased- 
infair (|Chiappini fc Matteuccil Il999f ). where in- 
fall of gas occurs at a faster rate in the inner- 
most regions than in the outermost ones ('inside- 
out' disk formation). Growing observational evi- 
dence for slow formation of disks in spiral galax- 
ies associated with continuing infall of primordial 
(or low m etallicity) gas over the lifetime of the 



disk (e.g. Braun fc Burton 1999f) seems to give 



suppo rt to the above scenario. iLubowich et al 
(|200fl[ ) demonstrated that continuous infall of low- 
metallicity gas is occurring in the Galactic Center. 

Of particular relevance to test these models 
are the abundances in the very outer galactic 
disks. Chemical evolution models for abundance 
gradients and the for mation of the Milky Way 



(jChiappini et al.l 120011 ) show that the steepness 



of the outer gradients are particularly sensitive 
to thresholds in star formation, to the halo-thick 
disk enrichment history, and to the radial variation 
of the disk formation timescales. They concluded 
that some abundance ratios increase substantially 
toward the outermost disk regions but that 'more 
observations at large galactocentric distances are 
needed to test these predictions'. The fact that 
N is almost consta nt with the galactocen tric dis- 
tance up to 18 kpc (Chiappini et aL 2003) refiects 
the high N production in AGB stars. The pri- 
mary N contribution from AGB stars is however 
very uncertain. At such large galactocentric dis- 
tances a N production threshold may be operating 
due to the paucity of late- type stars. 

Galactic chemical evolution also predicts that 
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the abundances of C, N, O, "c, and will 
be lower in clouds at the edge of the Galaxy 
than in any othe r clou d of the Galactic disk 
(jMaciel &: Quirezal Il999l ). The composition of 
these edge clouds should be similar to that of 
the early Galactic disk modified by infall from the 
halo. Thus the metallicity in clouds such as EC2 is 
expected to be similar to dwarf irregular galaxies, 
giving us an unique opportunity to study gas from 
the e arly stages of the formation o f the Galactic 
disk (jKobavashi fc Tokunagall200(J) . 

We have, therefore, searched for emission from 
a number of molecules towards EC2, concentrating 
on the northern CO peak (whic h is coinc ident with 
the T Tauri cluster identified bv lYasui"et al. 20061 , 
in order to constrain its physical conditions and 
chemical composition. Such molecular evolution 
in low metallicity Galactic clouds also provides a 
local laboratory in understanding molecular gas 
in extragalactic sources such as high-z quasars. In 
particular, since deuterated molecules are sensi- 
tive tracers of physical conditions and chemical 
pathways, we have also searched for a number 
of deuteriuni (D) b earing molecules, noting that 
Rogers et al.l (|2005l) have detected Di towards the 
Galactic anticenter. We also examine the role of 
SNR GSH 138-01-94 in providing the energy in- 
put for both the chemistry and star formation ob- 
served in EC2. 

2. Observations 

We have completed over 210 hours of obser- 
vations toward this object using the 12 m Ari- 
zona Radio Observatory (ARO) telescope (for- 
merly part of the National Radio Astronomy Ob- 
servatory), the 100 m Max-Planck-Institut fiir Ra- 
dioastronomie (MPIfR) Effelsberg telescope, the 
15 m James Clerk Maxwell Telescope (JCMT), 
and the 30 m Institut de Radio Astronomic Mil- 
limetrique (IRAM) telescope. The line observa- 
tions range in frequency from 4.83 GHz (H2CO 
li,o-li,i) to 492 GHz (Ci ^Pi-Spg). EC2 is ex- 
tended in CO and we have observed a number of 
different pos itions within ^40" of the three peaks 
seen in it bv lOigel et al. ( 1994h . Fig. [1] shows ob- 
served positions A through K on our contour map 
of CO 2-1 emission. Positions A, J and K corre- 
spond to po sitions copk l, copk 2 and copkS respec- 
tively from IPigel et al.l ( 1994 ) . In this work we 



present the results of our line search at position 
A (asooo = 02^48"" 38.5^ Jzooo = 58° 28' 28.1"), 
which at the time showed the strongest line spec- 
trum based on a twelve position HCO+ map. 
Our subsequent CO maps centered on position E, 
(02000 = 021^48™ 38. 5^ (52000 =58° 28' 58.3"), cov- 
ering the northern peak in emission, are also pre- 
sented. These maps revealed that the northern 
peak in CO emission was ~30" n orth-east of po- 
sition A, and ^25% stronger. See lLubowich et al 
(|2004l ) for observations of positions B and C 



Full details of observed line spectra for positions 
other than A, and CO maps covering the southern 
peak (centered at 02000 = 2*^ 48"' 23.6% J2000 = 
58° 23' 59.0") are available in lRuffld (|2006l ). 

Observations at the ARO 12 m were made be- 
tween June 2002 and January 2006 at 2 and 3 mm. 
Data were obtained in beam switching mode. The 
beam size (0b) of the 12 m telescope varies between 
93" at 72 GHz and 43" at 150 GHz. Observing 
conditions were generally good, with system tem- 
peratures ranging from 140 to 350 K for the 2002 
June observations. For subsequent observations, 
system temperatures ranged from 170 to 820 K. 
Dual channel SIS receivers were used in conjunc- 
tion with two spectrometer systems: Filter Banks 
1 and 2 (FBI and FB2), and a digital Millime- 
ter Auto Correlator (MACl). The 256 channel 
100 and 250 kHz filter bank spectrometers were 
used simultaneously, with 128 channels allocated 
to each polarisation, giving bandwidths of 12.8 
and 32 MHz respectively. The two polarisations 
were subsequently averaged during data reduction. 
The MACl was used to confirm filter bank data 
and set to a resolution of 50, 100, 200 or 400 kHz 
and bandwidths of 300 or 600 MHz. The ARO 
12 m produces spectral line data calibrated to the 
(K) corrected antenna temperature scale. 

We used the Effelsberg 100 m telescope in June 
2002, December 2002, January 2003, May 2004 
and November 2005 at 1.0, 1.3, 2, and 6 cm to 
observe transitions of SO, NH3, H2O, and H2CO 
with beamwidths (6'b) of about 30", 40", 55'" and 
160" and typical system temperatures of 400, 150, 
80, and 40 K on a main beam brightness temper- 
ature scale (Tmb)- Tinb to flux density conversion 
ratios were 1.6, 1.4, 1.8, and 2.4K/Jy, respec- 
tively. For the 1.0 cm SO lo-Oi line and the 2 cm 
H2CO 2i^i-2i^2 transition, single channel primary 
focus He'mT (High Electron Mobility Transistor) 
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receivers were employed. At 1.3 cm, for H2O and 
NH3 , a primary focus dual channel HEMT receiver 
was used, allowing us to add signals with orthog- 
onal linear polarization, thus reducing the effec- 
tive system temperature given above by ^^30%. 
At 6cm, H2CO li,o-li,i observations were taken 
with a four channel secondary focus receiver. Po- 
sition switching with appropriate offsets, alternat- 
ing between east and west in azimuth, was used for 
all measurements. An 'AK90' autocorrelator in- 
cluded four spectrometers with 2048 channels and 
bandwidths of 20 MHz each. Pointing was accu- 
rate to about 10" or better. 

In June 2004- July 2005, we used the JCMT to 
map the 2-1 and 3-2 transitions of several isotopes 
of CO and to observe the 492 GHz transition of 
Ci with beam sizes (^b) which range from 23" at 
220 GHz to 10" at 492 GHz. We also searched un- 
successfully for the HCO+ 3-2 transition at 267 
GHz. Data were obtained in a position-switching 
mode in a tracking coordinate frame with a chop 
throw of 600" at an angle of 270°. For the CO 
2-1 and HCO+ observations, the single channel 
A3 receiver was used yielding system temperatures 
from 260-330 K. For the CO 3-2 transitions the 
dual channel B3 receiver was used with two low- 
noise SIS mixers in single-sideband mode, yielding 
system temperatures from 420-800 K. Receiver W 
was used to observe the Ci transition in single side- 
band with a system temperature of 2400-4100 K. 
The JCMT produces spectral line data calibrated 
to the (K) antenna temperature scale, and this 
was converted to using the JCMT's rjiss (for- 
ward scattering and spillover) efficiencies. 

Finally, in October-November 2005 we used 
the MAMBO II bolometer at the IRAM 30 m 
telescope to make a 1.2 mm dust emission map 
of EC2. The bolometer comprises 117 pix- 
els with a resolution of 11". Zenith opacity 
ranged from 0.07-0.27 during the observations 
and the sky noise difference (between on and off 
source) was generally low, varying between 34 and 
90mJy/beam/2Hz. We made five maps of EC2 
centered on a2ooo=02:48:41.0, (52ooo=58:29:27.9, 
with an rms which ranged from 8.8-11.5mJy per 
beam. The five maps were combined using MOP- 
SIC. 



3. Analysis 

Detected lines, CO maps, dust maps and a 
summary of observations can be found in Figs. [T] 
through [7] and Table [TJ 

Before estimating abundances of the observed 
molecular species, we first have to determine the 
main physical parameters of the gas, i.e. its den- 
sity and temperature, taking into account the 
beam-filling by deconvolving the detected inten- 
sities. To estimate the latter we use the ^^CO 
and ^■^CO 2-1 maps in Fig. [5] and assume that 
the source, which has an elongated shell-like struc- 
ture, has a full- width at half power (fwhp) size of 
150" X 40" {0s = 77") for ^^CO, and 100" x 30" 
(6*3 = 55") for i^CO, C^^O and all other molecules. 
The peak intensities listed in Table[l]are corrected 
for these source sizes using T,„b ~ T^/r]i,{, with 
TTbf = ^s/(^s + ^b)i where T^bf is the beam filling 
correction factor. 

There are several tracers sensitive to both den- 
sity and temperature, though none that only 
traces density, but one that exclusively traces ki- 
netic temperature. The species allowing us to 
determine T^in is NH3. 

We have detected at least four of the five groups 
of hyperfine components in the (J, if) = (1,1) line 
of NH3 (see Fig. [21 the satellite features be- 
ing marked 'HF'). Line intensity ratios are con- 
sistent with optically thin emission. Gaussian 
fits of the groups of hyperfine components are 
(with increasing velocity) 23.3, 17.2, 80.2, 5.7 
and 11.1 mJy kms^^. For the (2,2) line we get 
38.0±8.0mJykms"^ The (2,2)/(l,l) line inten- 
sity ratio is then 0.28±0.07, with the error being 
dominated by the low signal-to-noise ratio of the 
(2,2) line. To derive total beam-averaged column 
densities in an inversion doublet, we use 



N{J,K) = 1.55 X 10^'' 



J(J+ 1) 



T^hdV (1) 



(e.g. eq. 1 in 'Henkel erall l200rf) that is valid in 
the optically thin limit. The column density N, 
the frequency and the integral are in units of 
cm~^, GHz, and K kms^^, respectively. As al- 
ready indicated in section [21 the conversion factor 
between the Jy and the T^b scale is 1.6. Thus we 
obtain total beam averaged column densities of 
2.52 X 10^2 and 5.21 x 10" cm'^ in the (1,1) and 
(2,2) inversion doublets. The (2,2)/(l,l) colmim 
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density ratio is then 0.21±0.05. With 

N{2, 2)/N{l, 1) = (5/3)e-^^/'=^ (2) 

and AE/k — 40.957 K, the difference in energy 
between the (f ,f ) and (2,2) levels, we derive a ro- 
tational temperature T = 20 ± 3 K. 

With the also populated inversion singlet (0,0) 
leveQ at 23.19K below the (l,f) state, the total 
NH3 column density, summed over the (0,0), (l,f ) 
and (2,2) levels is A^(NH3) = 5.77 x lO^^cm'^, 
with an uncertainty of ±f 5 per cent, averaged over 
the 40" beam. J > K levels require e xtremely 
high densities (jMauersberger et al.l Il985^ and are 
therefore not relevant. 

Using this rotational temperatur e, which should 
be close to the kinet i c temp erature ()Walmslev fc Ungerech 

the density of EC2 



160/I8O = 500 (|Frerking et al.lll982f) . 

For the rotational HCN 1-0 transition, our de- 
tected hyperfine structureQ is compatible with op- 
tically thin emission, implying that other lines 
should also be optically thin (HCN 1-0 is among 
the optically thickest 3 mm lines after CO). There- 
fore, we have computed LVG models for several 
other species, and start with those for which we 
have measured more than one transition. For 



1983t iDanbv et all Il988l) 



SO, using the coUisional rates of iGreenI (|1994| ). 
we can fit the Iq^Oi and 32-2i lines simultane- 
ously, but only with n(H2) = 1-4 x 10^ cm""^, 
an order of magnitude larger than that derived 
from H2CO and two orders of magnitude higher 
than that obtained from ^^CO, implying the pres- 
.£iice of small-scale clumping. The CS 2-1 and 3- 
W deconvolved line temperatures ca n be roughly 



was estimated using Large Velocity Gradient 
(LVG) models for a number of molecular species, 
for a spherical cloud with constant density (see 
iHenkel et ai]ll980l : [Maiiersberger et al.lll990L and 
references therein) . The results from these models 
are included in Tabled First, we use a LVG model 
involving 40 ortho-H2CO stat es up to K„ = 3 and 



the collision cross sections of iGreenI (|1991l ). We 
find that the observed deconvolved line intensities 
for both the cm and mm lines are well reproduced 
with a density n(H2) = 1.2 x lO'' cm~'^. Second we 
fit the "CO and C^^O spe ctra to a n LVG model 
using coUisional rates from iFlowed ( 2001 ) and an 
adopted ortho-para ratio of 3 for H2 , although the 
results are only weakly dependent on this. We first 
fitted the three ^^CO deconvolved lines to derive a 
density of n(H.2) = 3.2 x 10^ cm~^. With this den- 
sity, C^*0 (only two lines) and ^^CO could then 
also be fitted. This gives isotopical abundances of 
approximately 100:10:1 for i2C0:i3C0:Ci80 (see 
below). This can not reflect the actual isotopic 
abundance ratios, as the optically thin ^^C/^^C 
ratio is ~89 in the solar system and 60 t o 70 in 
local diffuse clouds ( Wilson fc Roodl[l993 ). Using 
nivg(X)/ri(H2) from Table [21 the best fit also gives 
a fractional abundance with regard to H2 of ~10~^ 
for I^CO, weh below the value of 2 x 10"*^ typi- 
cally found in local metal-rich molecular clouds. 
For the C^^0/H2 ratio the model gives -10~^ 
compared to the local value of 1.7 x 10~^ with 



fitted (using coUisional rates from [Turner et al 



19921 ) with the density of 1.2 x lO^'cm-^ obtained 
for H2CO. The model reproduces HCO+ 1-0 de- 
convolved line temperatures for both the H2CO 
and the SO densities (using coUisional rates from 
Flower Il999f) . The HCO+ model line intensities 
appear to be consistent with the observed J — 3-2 
upper limit, at least in the case of the lower den- 
sity of 1.2 X 10"* cm"'^. For HCN 1-0 we summed 
the deconvolved integrated intensities of the three 
hyperfine components and then divided by the 
average of the line widths to get a peak inten- 
sity. Using coUisional rates from Schoie r et al 



(|2005l ). this peak deconvolved line intensity was 
fitted for both the H2CO and the SO densities 
(1.2 X 10** cm^^ and 1.4 x 10^ cm~^ respectively). 
The higher density model is likely to be the cor- 
rect one; as in the case of HCN, radio observa- 
tions are mainly sensitive to high density gas. The 
model reproduces HNC deconvolved line temper- 
atures for both the H2CO and SO densities (using 
rates from lSchoier et a"Lll2005l ). N2H+ was not de- 
tected, which is highly unusual and presumably a 
sign of N depletion. 

There are several ways to explain the discrep- 
ancies between the expected and derived from 
observations ^^CO:^'^CO:C^^O abundance ratios, 
and the surprisingly low fractional abundances of 
i^CO and C^^O: (i) optical thickness of ^^CO; (ii) 



-'The (0,0) level is ortho-ammonia, with a statistical weight 
of 2, but as there is only one state, this is factored by 0.5. 



^The HCN spectra in Fig.[3]shows spikes near — llOkms"^, 
increasing the strength of what should be the weakest hy- 
perfine component. Spectra taken 30" from position A in 
2005 Dec. confirm the hyperfine structure for HCN. 
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^^CO enhancement by fractionation in cloud edges 
due to charge exch ange reactions with C"*" (e.g. 



Watson et al.l Il976f l: or (iii) small scale clump- 
ing, suggested by the H2 density derived from 
the SO line intensity ratio. The energy barrier 
for the fractionation r eaction ^^C"*" + ^^C O < — > 
is ~35K (jWatson et al.lll97(ih . which 
requires low temperatures for efficient fractiona- 
tion, although at 20 K the enhancement would 
not be large. Low Ay and a high UV field are 
also required for the production of C"*". However, 
for higher UV fields, since ^^CO self-shields more 
effectively than ^"^CO, selective photodissociation 
destroys ^^CO faster than ^^CO. So depending on 
the respective rates of the fractionation reaction 
and photodissociation, ^^CO fractionation may or 
may not be taking place. Given the range of 
densities derived, small-scale structure and opti- 
cal depths effects are the most likely explanations 
for the discrepancies (section [S] examines the effect 
of high UV fields and low Ay in greater detail). 

In order to derive accurate molecular abun- 
dances one needs multiple transitions of a vari- 
ety of isotopes to account for optical depth and 
excitation effects. Furthermore, the use of opti- 
cally thin transitions of isotopes to derive accu- 
rate abundances of the main species is also fraught 
with difficulty since the underlying isotopic ratios 
may be different in this gas at the edge of the 
Galaxy than in the more local molecular clouds. 
To derive molecular abundances for species other 
than CO, H2CO and NH3, we have assumed LTE 
and an excitation temperature of 20 K. Most of 
the derived abundances are not very sensitive to 
temperature, typically varying by less than ±10 
per cent for T = 20 ± 3K, except for the abun- 
dances derived for CO 3-2 (±22%), Ci (±18%), 
SO (±13%), H2CO (±12%) and NH3(±12%). 

For optically thin line emission from linear 
molecules, we find that the total column density, 
N (cm~^), of a molecule with normalised electric 
dipole moment p, (in units of 10"^* esu) with an 
excitation temperature Tox emitting at frequency 
v in a, transition J— > J-1 can be written in the form: 



detections we use the 3cr upper limits for the inte- 
grated intensity: 



N 



1.603 X 10^3 T, 



<=''g(J+l)To/2T„ 



rp2 
^0 



r„,bdu. 



where Tg = hv/k and the deconvolved inte- 
grated intensity for a Gaussian line profile is given 
by 1.06TnibAy in units of Kkms^^. For non- 



A(Kkms"^) = 3rTOs,„b(l + a)VAz^AU, (4) 

where a is the calibration uncertainty taken to be 
0.15, Ai/ and AV are the spectral resolution and 
expected line-width (2kms~^ is adopted here), 
both measured in kms~^. For optically thin emis- 
sion, the fractional abundance varies proportion- 
ally to the integrated line intensity. For CN we 
summed the intensities of our two observed hyper- 
fine components and adopted a relative intensity 
of 0.456 to determine the integrated intensity and, 
thus, N{CN). For HCN the sum of the three hy- 
perfine components was used, while for C2H we 
used the strongest hyperfine component at 87.317 
GHz and adopted a relative intensity of 0.4167. 

We take a slightly different approach for SO, 
which has a '^S ground state. We cal culate A^(SO) 
directly using the formula given bv iTakano et al 



(|l995l) with the partition function from the JPL 



Mole cular Spectroscopy Catalogue (jPickett et al 
1998[). For Cl w e use a similar approach to that 



Tauber et all (|l995f) . To estimate iV(H2C0) 



of 

we use the LVG modelled density: iV(H2C0) = 
n(H2C0) X Ipc X AU, where ^(HaCO) is the 
best fit abundance to the LVG calculation (1.4 x 
10-6 cm-3) and AV = 1.3 kms^^ (the li^o - 
line is widened by hyperfine structure) . Table [1] 
lists the derived column densities. 

Table [2] compares line temperatures and column 
densities for EC2 derived from our LVG models 
and observed deconvolved integrated intensities, 
and it can be seen that in most cases there is good 
agreement between N„-±,{X) and A^ivg(X). The 
LVG column densities were calculated from the 
best fit model densities with iVivg(X) — nivg(X) 
x Ipc x AV. It should be pointed out that our 
assumption that Tex = Tkin = 20 K could be con- 
sidered unrealistic, since excitation temperatures 
are crucial for abundance determination. CO is 
almost certainly thermalised, but other molecu- 
lar species may not be thermalised with similar 
certainty. At the low densities we have derived 
from the observations, molecules like HCO+, HCN 
and HNC would not be thermalised. However, 
our LVG results provide a realistic estimate of the 
real excitation temperature (Tex in Table [5]) for 
a range of our detected molecules. Recalculating 
with these specific values (instead of Tox — 20 K), 
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gives column densities that only differ by factors 
of between 0.5 and 3.7 (except CS 3-2 by 8.2). 

To aid comparison with local molecular clouds, 
we have derived molecular abundances relative to 
that of HCO"*", assuming that the 1-0 transition is 
optically thin. Our upper limit to H^'^CO+ gives 
a lower limit to the ^^C/^'^C ratio of around 17, 
much less than t he lower limit of 201zfcl5 der ived 
from CO 1-0 by IWouterloot fc Brandl (|l996l ) to- 
wards WB 89-437 at a galactocentric distance of 
16.4 kpc, although we note that the value they 
derive from the 2-1 transition is much smaller, 
104±60. Table [3] shows the result and also the 
equivalent abundance ratio for the nearby dark 
clouds L134N (center position), TMC-1 (aver- 
age), and the range of values derived in L134N 
(|Dickens et al.l l2000t IPratap et all [Im^). Note 
that these authors use isotopic species to deter- 
mine the column densities of optically thick lines. 
In particular, they adopt ^'^C/^^C = 64 for deter- 
mining the column densities of HCO"'", HCN and 
HNC and ^^O/^^O = 500 to determine the CO 
and "'^'^CO column densities. Finally we note that 
they used a dipole moment of 4.07 D for HCO+ 
compared to the value of 3.9 D used here. Use 
of our value would increase their abundance ra- 
tios by about 10%. Table [3] also gives fractional 
abundances relative to molecular hydrogen using 
-^(H2) = 7.4 X 10^^ cm~^ and compares these val- 
ues with_Jhosefor the translucent clouds observed 
bv lTurneil (|2000l ). 

Dust emission can be used to estimate the mass 
of cold molecular clouds and we use our detected 
peak continuum emission at 1.2 mm to calculate 
the opacity and dust-to-gas ratio in EC2. First, 
we estimate Tj^, the dust opacity at frequency v, 
using 

T, = B(iy)/7r(0b/2)'B,(Tdust), (r « 1) (5) 



relation pdust — t^/^vL^ and it therefore follows 
that the dust to gas ratio is: 



from iHildebrandl (|l983h . where B^iT^ust) 

is the 

Planck function, Tdust = 20 K, the observed dust 
emission B{i') = 20.3 x 10"'^ Jy, and the beam size 
6'b = 11", which gives a value of = 3.25 x 10^**. 
The dust opacity per unit mass column density at 
frequency v, is given by — SQ^/iap, where Q^, 
is the dust emissivity which varies as i^^ [5('^) 

B^{Tdust)], the average dust grain radius a = 
0.1 /xm, and the density of the dust p = 3g cm~'^. 
Value s for Q,^, a and p a re tak en from Hilde brandl 
(|l983l) . From IWhittetl (|l992l . pl50) we note the 



K,yL n(H2)^mH ' 



(6) 



where L = DOh is the size of the source and D 
is the heliocentric distance (14-20 kpc). Tabled 
lists derived dust to gas ratios for D — 14 and 
20 kpc and a range of values for /3 from 1 to 2 
(taken to characterise dust grains that range from 
amorphous to crystalline), as well as appropriate 
values for (3i2oo and k^. 

It now remains to relate the dust to gas ratio 
to possible values of extinction in EC2. For the 
standard dust-to-gas mass ratio of 0.01 in the local 
ISM, it is found that NCR)/E(B-y) = 5.8 x lO^i 
atoms cm^'^ mag~^ ( Bohlin et al. 1978h and with 
Ry = Ay/E{B — V), we can write: 



5.8 X 1021 



gas 



0.01 



(7) 



As discussed in iRuffle et al.l (|2004h . the value of 
Ry may be lower in certain regions of the ISM, 
so we calculate values of Ay for Ry — 3.1 and 
Ry — 2.0 at the position of maximum contin- 
uum emission. As our CO maps show, EC2 is 
clumpy in structure and the peak dust detection 
of 20.3mJy/beam (0b = H") traces the peak CO 
emission of one such clump. Over the map, the 
rms value of 6.35mJy/beam is ~3 times less than 
the peak, so extinction in EC2 as a whole is likely 
to be correspondingly lower. 

Finally, we estimate upper and lower limits for 
the mass of EC2 based on a density of n(Il2) = 
1.2 X lO^'cm"^. From our maps of the two CO 
2-1 peaks in EC2 (Fig. [1]) we estimate that the 
cloud consists of two dense clumps ~75" in di- 
ameter, which yield a mass range of 9.3 x lO'' 
to 2.7 X 10^ Af0 for the distance range 14-20 kpc. 
Taking the clumpy structure of EC2 into account, 
the average density may be of the order n(H2) = 
10'^ cm~'^, so our mass estimate is Mf,c;2 ^ l?^"^? 
(based on their CO detections IPigel et al.l ( 19941 ) 
estimated Mec2 = 3.7 x IO^Mq). This value 
can be compared with a virial mass estimate us- 
ing the method of iMacLaren et al] (1988, equa- 
tion 3), Mvir/Mo = k2rAV'^, where = 190 
(for p cx l/r), r is the radius of the clumps in 
pc, and AV is the ^^CO line width (FWHM) in 
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kms~^. iMacLaren et alj argue that the virial the- 
orem will only yield a meaningful result if the esti- 
mate of AT/ is found from the average throughout 
the whole volume of the cloud, including central 
core regions. An optically thick line such as ^^CO 
will only probe the cloud to the surface where the 
optical depth is unity. On the other hand, an opti- 
cally thin line, such as ^'^CO should allow sampling 
of velocities throughout the central, denser regions 
of the cloud. However, taking the "radius" of EC2 
from the ^^CO observations should avoid the ef- 
fect of excluding the lower density outer regions, 
i.e. an overestimation of the total mass. For the 
distance range 14-20 kpc, the radius of each of the 
two clumps, r is 2.5-3.7pc and the line width for 
the strongest ^^CO 2-1 transition is 1.87 kms~^, 
which yields a total virial mass Mvir = 3.3-4.9 
X 10^ Mq, somewhat less than our direct estimate 
above. The lower virial mass suggests that the 
cloud is gravitationally unstab le and a likely site of 
star formation as suggested bv lYasui et al. [ ((20061 

4. Discussion 

The data in Table [3] indicate some interesting 
differences between molecular abundances relative 
to HCO+ in EC2 and local dark molecular clouds. 
In particular sulfur-bearing molecules, (CS, SO), 
appear to be very over-abundant; the nitrogen- 
bearing species HCN is comparable while HNC 
appears to be marginally under-abundant - the 
HNC/HCN ratio is less than one, whereas in dark 
clouds it is usually close to or slightly greater 
than one - and the radicals CN and C2H are very 
much over-abundant. The latter result is typi- 
cal of photon-dominated regions in which photo- 
processes ensure that radicals exist in high abun- 
dance. At first glance it would appear that the 
molecular clouds at the edge of the Galaxy would 
be less likely to show the effects of PDR chemistry; 
certainly the flux of UV photons must be much 
less at the Galactic edge than in the local ISM. 
However, the critical parameter for determining 
whether or not photons dominate chemistry is the 
ratio of UV flux to grain surface area, particularly 
the area of the small grains which extinguish UV 
photons most efficiently. At large galactocentric 
radii, the metal abundances relative to hydrogen 
are expected to be much reduced, consistent with a 
reduced dust-to-gas ratio, as discussed previously. 
In addition, although the region of EC2 appears 



to contain young T Tauri stars ( Yasui et al.ll2006f ). 
there is no evidence of the late-type stars which 
produce dust grains. 

5. Chemical Modelling 

To investigate the properties of EC2 and to 
see whether it is typical of material which has 
been less processed, we have made a pseudo-time- 
dependent chemical kinetic model. This uses 
observationally derived temperatures and den- 
sities and varies elemental abundances, photon 
flux, cosmic ray ionisation (CRI) rate and gas- 
to-dust ratio in an attempt to fit the observed 
results. Chemical reaction rates were taken from 
Rate99: The UMIST database for astrochemistry 
()Le Teuff et al.ll2000l) . 

Unless otherwise stated, the model assumes 
Tkin = 20 K, n(H2) = 1.2 x 10'* cm'^, based on 
NH3, H2CO, and GO observations, and uses the 
initial elemental abundances listed in Table O For 
the metals (G, N, O, S, Si and Fe) these are re- 
duced by a factor of five from those typically used 
for local ISM clouds; we find that without this re- 
duction in elemental abundances, agreement with 
observations is worsened for Ay > 1 mag. We have 
investigated varying the GRI rate and the visual 
extinction (Ay) in order to try and find the best 
fit to the observations. We also varied the UV 
photon field, initial abundances, and, specifically, 
the abundances of sulfur or nitrogen. 

Our JGMT GO maps probed at higher resolu- 
tion (15-23"), but the bulk of our detections (with 
the ARO 12 m) covered beams of 42-70". Our Ef- 
felsberg detections had beams of 30" for SO, 40" 
for NH3, and 55" and 160" for H2CO. As described 
in section [31 our detected peak intensities were de- 
convolved for these beam sizes (^b) and our esti- 
mated fwhp clump size {9s)- The peak GO 2-1 
emission seen in our maps is within ~30", or half 
a typical beamwidth, of position A. Our ARO 12m 
detections of GO, CS, SO, HGO+, HCN, C2H and 
C 3H2 at nearb y positions B, C, D, E, G, H and 
J ()Ruffldl2006l) are commensurate with their posi- 
tion A detections, given that the intensity contours 
seen in Fig. [1] would be 'smeared' by the relevant 
beam sizes (55-70"). At a distance of 20 kpc, a 60" 
beam subtends a linear distance of about 6 pc or a 
column density of about 2 x 10^'^ cm~^ for n(H2) = 
1.2 X 10^ cm~'^. Since our GO maps show that EC2 
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is extended and clumpy, the 'one-point' model 
adopted is simplistic; nevertheless, given the lack 
of data on very small scales, one might hope that 
the model helps constrain some of the global pa- 
rameters in the cloud in a 'beam-averaged' sense. 
321 models with varying parameters were run and 
the results of testing each one agai nst our obser- 
vations are detailed in Ruffle ( 20061 ). 

To test the agreement between each model and 
our observations, we looked at the ratio of the ob- 
served abundance relative to H2 from Table [3] to 
the model prediction (or the inverse of that ratio 
if it was greater than 1) for that species. For each 
molecule/transition we then square its agreement 
factor and sum these squares for each model, nor- 
malising to get a 'Fit' number between and 1. 
This should be closest to 1 for the model which 
agrees best. In the case of our multiple transition 
detections of ^^CO, ^^CO and C^^O, we have used 
a a weighted mean of each isotopomer's observed 
column densities, where a — T^b/rmSmb- How- 
ever, as the models only consider the main isotope 
of an atom (with the exception of H and D) we 
have assumed that the isotopic ratios have their lo- 
cal interstellar values (^^C/^C = 60, G^^O/C^^O 
= 500) and factored our observed (mean) abun- 
dances relative to N{}i2) accordingly. For the 
"'^^CO and C^^O transitions, this gives us a rea- 
sonably consistent value of N{C0)/N(R2) ~ 1-4 x 
10"^, although the main CO lines give a value 
closer to 10"^. Since the main lines of ^^CO are 
optically thick (0.6-2.2 from Table [2]), the ratios 
derived from observation are lower limits. 

We have selected 10 models for detailed exam- 
ination, which can be compared in Tables El [71 [8] 
and O We investigated the agreement factor at 
steady state which for our choice of parameters is 
reached, in general, at times less than than 10^ 
yr. Fig. |8l plots 3D parameter surfaces of these 
test results for various combinations of the models, 
where two parameters are varied and the others 
are kept constant. For ease of comparison, num- 
bered disks, marked on the 3D plots, indicate the 
examined models relative to the parameter sur- 
faces. Fig. [9] shows fractional abundances varying 
over time for five of the models. 

Varying CRI and Ay 

We investigated the agreement tests for values 
of Ay between 0.5 and 10 mag and a CRI rate 



between 0.5 and 30 times the standard interstellar 
rate of 1.3 x lO'^'^s'^ The UV field was set to 
the local ISM value. Although many of the mod- 
els give very similar results, the agreement fac- 
tor is highest for a CRI rate above 10 x the ISM 
value and an Ay of 1 mag. Fig. [8^ shows the gen- 
eral trend of varying Ay and CRI. It can be seen 
that increasing the CRI rate beyond 20 x the ISM 
value extends agreement in the model to Ay = 
2 mag. What is striking however, is the failure of 
the model for values of Ay above 2 mag. It should 
be noted that model results for Ay < 0.5 mag 
should be treated with caution, as the model may 
fail for such low values of Ay, since it does not 
treat the self-shielding of II2 and CO accurately. 

Table [7] shows how well these models actually 
agree with the observations for four selected mod- 
els (1-4). Model 1 uses 'standard' ISM values for 
the Ay and CRI rate, Model 2 is the best fit model 
(CRI = 2.6 X 10-16 s~\ Ay = 1 mag), while 
Model 3 uses CRI = 1.3 x lO^^^ s^^ and Ay = 
2 mag, to illustrate some effects which arise from 
changing these parameters. Model 4 is the same 
as Model 3, but with a UV field 10 times the local 
rate. 

For Model 1 the predicted CO/H2 ratio is 1 or- 
der of magnitude higher than all the observations. 
Models 2 and 4 agree within a factor of 2 for ^^CO 
but underestimate ^^CO and C^^O by ~1 order of 
magnitude, whereas Model 3 reverses those agree- 
ments. Model 1 does badly at reproducing the 
observed value for atomic C, since when the Ay 
is 10 mag the cloud is shielded from UV radiation 
and the bulk of the carbon becomes incorporated 
into CO. When the Ay is lowered and/or the CRI 
rate increased, however, CO can be broken down, 
so that the abundance of Ci remains significant 
even at steady-state. We find that an Ay of 1 
or 2 mag gives a better agreement with the EC2 
observation of N{C)/N(Ji2)- 

The predicted value of C2II/II2 in Model 1 
agrees with the observations to within a factor 
of 2, while the other three Models are 10 to 40 
times too high. Conversely, Models 1 and 3 over- 
predict the CN/H2 ratio by over 2 orders of mag- 
nitude, with Models 2 and 4 over-estimating it by 
40 and 50 times respectively. The time-dependent 
results show that Models 2 and 4 have good agree- 
ment overall once t is greater than about 2000 yr, 
whereas neither Model 1 or 3 show good agreement 
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beyond the first 500 years (see Fig. [9f). 

Varying UV photon field, CRI and Ay 

We have also looked at the effect of varying the 
photon field in the models. Fig. [H}3,c shows agree- 
ment factors for models with varying Ay and CRI 
for (b) a UV field increased by a factor of 10 com- 
pared to the local ISM value, and (c) a UV field 
reduced by a factor of 10. It can be seen that a 
higher UV field improves agreement up to an Ay 
of 2 mag regardless of CRI rate (although for CRI 
20 X this is extended to 3 mag), with Model 4 
the best fit. A lower UV field produces very poor 
agreement (e.g. Model 5). 

Fig. [8]i plots fit results from models where the 
CRI rate is fixed at 20 times the standard value, 
Ay is varied, and the UV field is varied from one- 
tenth to 80 times the local interstellar field. As 
expected, an increase in the UV photon field ex- 
tends agreement to higher values of Ay (4 mag for 
80 X UV). Table [8] shows the agreement for three 
of these models. Interestingly, 20 x UV produces 
a good fit for atomic carbon at an Ay of 1 mag 
(Model 6), whereas 40 x UV at an Ay of 3 mag 
gives a good fit for CO (Model 7). At steady-state 
CO is primarily destroyed by photons (only if the 
UV field is strong enough, or the Ay low enough) 
or by He"*" and H|^, whose abundances increase 
with increasing CRI rate. Again a UV field below 
the local value produces very poor general agree- 
ment, but apart from very high UV fields (~70 x) 
at an Ay of ~4 mag, there is poor agreement for 
any Ay above 3 mag. Interestingly the high UV 
(20 X ) Model 6 shows remarkable agreement from 
the earliest times and Model 7 (40 x UV) agrees 
well from t = 2 x 10** yr. 

The observed upper limit to the DCO+ /HCO+ 
ratio, 0.08, is not a useful limit since all models 
with the exception of Model 1, satisfy this con- 
straint. For C2II/II2 and CN/II2, the ratios are 
overestimated by 2 orders of magnitude in the low 
UV Model 5, but the fit improves with increasing 
UV field, with good agreement reached at 20 x 
UV (Model 6). Finally, we note that there is a 
much better fit to the ammonia abundance in low 
UV Model 5 (which is also achieved for NH3 in 
Model 3 with UV = 1 and CRI = 10 times the 
local ISM rates). 



Varying initial abundances and Ay 

Initial elemental abundances (see Table [5]) were 
varied between local ISM values and 10 per cent 
of those values. The parameter surface in Fig. [5^ 
shows that for Ay above 1 mag, the agreement 
worsens appreciably for initial abundances (lA) 
above 20 per cent of local cloud values. Table [9] 
compares Models 8-10, and for lA = 0.1 (Model 8) 
there is good agreement, except for atomic C and 
CN (~1 order of magnitude too great), NH3 (^1 
order to little), and SO (underestimated by 5 or- 
ders). Model 8's agreement varies considerably in 
the first 10,000 years or so, before setthng down 
at a steady-state below Model 2 (lA ~ 0.2). 

The observed CS(2-1)/H2 ratio is fairly well 
matched by our best fit Models 2 and 4 and very 
well matched by Models 7 and 8 (UV = 40 and lA 
= 0.1 respectively). Except for Model 1, though, 
the predictions for SO are too low by several or- 
ders of magnitude. Therefore, we tested the effect 
of increasing the initial abundance of sulfur rela- 
tive to other initial abundances (which are set at 
20 per cent of local cloud values) . As expected the 
abundances of sulfur-bearing molecules scale lin- 
early with any change in the initial abundance of 
S. For the N-bearing species, all models except 1 
give good agreement for HC3N and N2II+, but the 
derived abundances are upper limits, so this result 
must be treated with some caution. For HCN and 
HNC the picture is less clear, but abundances are 
again well matched for our best fit Models 2 and 
4. Models 3 and 5 are in good agreement with 
the NII3/H2 ratio, while all other models (except 
Model 1) predict ratios 2 to 5 orders of magnitude 
too low. For CN the high UV Model 6 agrees with 
the observed abundance. 

Varying n(Il2) and Ay 

Fig. [5f shows the effect of varying n(Il2) and 
Ay with a UV field set to the local ISM value and 
a CRI rate 20 times the local ISM value. Densities 
above n(H2) = 1.2 x lO^'cm^'^ produce a marked 
fall-off in agreement. We also found that, even for 
a UV field 10 times the local interstellar value, the 
visual extinction must be below Ay = 3 mag for 
reasonable agreement. 
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Deuteration 

We failed to detect any deuterated species 
in this survey, with the upper hmits being: 
DCO+/HCO+ < 0.08; DCN/HCN < 1.5; C2D/C2H 
< 0.06 and HDCO/H2CO < 0.01. The models pre- 
sented here use the underlying D /H ratio from the 
local ISM. Model 1 predicts DCO+/HCO+ ~0.25, 
which is higher than that seen in local dark clouds 
such as TMC-1, a result of the lower abundance of 
heavy elements. Models 2 and 4, however, predict 
very low DCO+/HCO+ ratios (--lO^^)^ whereas 
Model 3 predicts ~10^^. The large upper limits 
on DCN/HCN and C2D/C2H mean that these are 
not useful in constraining the models. Likewise for 
HDCO and DCO"*" the fractionation ratios agree 
except for Model 1. 

Fractional abundances over time 

To this point our model results have been dis- 
cussed in terms of steady-state abundances and 
it is pertinent to ask whether there might be 
any significant difference if we consider the time- 
dependent results. Fig. [§] plots the fractional 
abundances over time for Models 2, 4, 6, 7 and 
8. For our best-fit model, Model 2, steady-state 
is reached very quickly, after around 5,000 years 
(apart from some slight changes in HCN and 
HCO"*"). Depending on the choice of the UV field 
and visual extinction, other models reach steady 
state on somewhat shorter (~500 yr. Model 6) or 
somewhat longer (~ 10^ yr. Model 7) times. The 
longest time to reach steady state is ~5 10^ yr 
(Model 8) which has the lowest initial elemental 
abundances. Given that the likely age of the EC2 
molecular cloud is greater than a few thousand 
years, we expect that the steady-state results can 
be used to compare with the observed abundances. 

6. Conclusions 

In order to deduce its physical and chemical 
properties, we have observed continuum emission 
and a large number of molecular transitions in 
EC2 at the Galactic edge. We have also made 
CO maps of EC2 and used these to calculate de- 
convolved line intensities. A temperature of 20 K 
was estimated from hyperfine detections of am- 
monia and a gas density of n(H2) ~ 10^ cm^'^ 
was determined by comparing LVG models of a 
number of species to their deconvolved line de- 



tections. Molecular abundances were also deter- 
mined from the LVG models and found to be in 
good agreement with abundances calculated di- 
rectly from the deconvolved line intensities. From 
the peak continuum emission we calculated a dust 
mass for EC2 > 150 M0 and a dust-to-gas mass 
ratio > 0.001. 

Through the use of chemical models we have 
been able to establish the most likely chemical and 
physical properties of EC2 (although we note that 
not all observed abundances can be reproduced in 
a self-consistent manner - see Table |6] where our 
chemical model agreement results are ranked by 
fit). There is an indication that heavy elements 
may be depleted by about a factor of five rela- 
tive to local molecular clouds (similar to those in 
dwarf irregular galaxies and damped Lyman al- 
pha systems). Such reduced abundances may be 
attributed to the low level of star formation in 
this region and are probably also related to the 
continuing infall of primordial (or low metallicity) 
halo gas since the Milky Way formed. The models 
also suggest a high UV photon field in EC2 (10- 
20 X local values), where an increased UV field 
allows for values of up to 4 mag, especially if 
this increased field is combined with an increase 
in CRI (10-20 x), although the models are less 
sensitive to increases in the CRI rate. High CRI 
rates (>20 x the ISM value) without an increase 
in UV field only allow for extinction up to 2 mag. 
Gas densities much above n(H2) = 1.2 x 10'*cm~'^ 
are excluded by the models, even if the UV field 
is increased. Some of our models indicate that 
steady-state is reached very quickly after around 
5,000 years and that a high UV field can reduce 
this time to just ~500 yr. For a very high UV pho- 
ton field (40 X ) and an extinction of Ay — 3 mag, 
however, steady-state is not reached until 10'' yr. 

In the context of ratios relative to HCO+, 
sulfur-bearing molecules appear to be very over- 
abundant by at least an order of magnitude com- 
pared to local dark clouds. The observed high 
abundances (again relative to HCO+) of the radi- 
cals C2H and CN are typical of photon-dominated 
regions (PDRs). This may be related to a large 
value of the UV flux to grain surface area when 
compared to local clouds. In particular, we find 
that our best-fit models are consistent with re- 
duced elemental abundances and a low dust-to- 
gas mass ratio. In addition, although EC2 does 
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contain young stars, there is no evidence of the 
late-type stars which produce dust grains, thereby 
justifying the assumption of a high ratio of UV 
flux to grain surface area. We conclude therefore, 
that despite the position of EC2 in the Galaxy, 
UV photons (rather than cosmic rays) play an im- 
portant role in establishing its detailed chemical 
composition. The observed clumpy structure also 
leads to an enhanced role for PDR-like chemistry 
in EC2. 

Given that EC2 is in a region of extremely low 
gas pressure and very small spiral arm pertur- 
bation, the question remains as to the origin of 



the structure and chemistry in EC2. IStil fc Irwin 
i2001) speculate that old SN shells may be a 
source of dense clouds in low density environments 
such as the outer Galaxy. They go on to show that 
the SNR associated with EC2, GSH 138-01-94, is 
the largest and oldest SNR known to exist in the 
Milky Way. It consists of a H i shell at a kinematic 
galactocentric distance of 23.6 kpc, with an expan- 
sion velocity of 11.8±0.9kms^^, an expansion age 
of 4.3 Myr and a timcscalc for dissolving into the 
ISM of 18 Myr. ^til fc Irwini associate EC2 with 
the approaching side of the Hi shell, reducing the 
distance of EC2 to i? ~ 23.6 kpc, compared with 



28 kpc ([Digel et al.lll994[ ) , and in closer agreement 
with the photometric distance range of 15 -19 kpc 



for the B star MRl (jSmartt et al.ill996f ). EC2 



could even have been formed by GSH 138-01-94 
from swept-up interstellar gas, through Rayleigh- 
Taylor instabilities. It could, therefore, be as 
young as the ages derived from the time-dependent 
calculations, discussed above. Although shock 
chemistry may be driven by the SNR, the short 
time-scales to reach steady-state imply that UV- 
driven chemistry can reset shock abundances. One 
test of the importance of shock chemistry would 
be to search for SiO emission, a result of an en- 
hanced Si gas-phase abundance due to sputtering 
of interstellar grains. 

We conclude that the formation, structure and 
subsequent chemistry of EC2 may be the direct 
result of shock fronts from GSH 138-01-94 prop- 
agating through the medium sometime between 
1,000 and 10,000 years ago. The role of GSH 
138-01-94 i n rela tion to EC2 is reinforced by 
Yasui et al. I (|2006l) . who conclude that the three- 
dimensional geometry of GSH 138-01-94, EC2 
and the embedded star cluster, as well as the clus- 



ter's age ('--^l Myr, much less than the SNR's ex- 
pansion age of 4.3 Myr), strongly suggest that that 
the cloud collapse and subsequent star formation 
observed in EC2 was triggered by the SNR Hi 
shell. 
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Table 1 

Summary of observations and detections toward EC2 position A. Peak intensities 

CORRECTED FOR SOURCE SIZE USING Tmb = T^hhl, WHERE rybf = 9^ / {9l + O^)] COLUMN DENSITIES AND 
abundances estimated FOR AN EXCITATION TEMPERATURE OF 20 K (SEE SECTION . 



Molecule 


Transition 


Freq. 


Tmb 


{km s""*") 


l\V 

(kms~ ) 


rms 
(K) 


(cm"2) 




1-0 


115.271 


7.863 


2.70 


0.650 


0.151 


2.48 ± 0.06 X 10^'^ 


"CO 


1-0 


110.201 


1.499 


2.24 


0.272 


0.131 


4.24 ± 0.45 X lO^'' 




1-0 


109.782 


0.119 


2.05 


0.683 


0.027 


3.11 ±0.94 X 10" 


c"o 


1-0 


112.359 






0.667 


0.012 


< 5.62 X 10" 


CO 


2-1 


230.538 


6.289 


2.76 


0.406 


0.026 


8.81 ±0.05 X IQis 


i^CO 


2-1 


220.399 


1.897 


1.87 


0.425 


0.073 


1.90 ±0.09 X 10^^ 


C^SO 


2-1 


219.560 


0.167 


1.48 


0.427 


0.036 


1.33 ±0.41 X 10^"' 


CO 


3-2 


345.796 


3.715 


2.56 


0.542 


0.020 


4.92 ±0.03 X lO^s 


i^CO 


3-2 


330.588 


0.716 


1.69 


0.567 


0.027 


6.37 ±0.33 X lOi"* 


CiSO 


3-2 


329.330 






0.569 


0.022 


< 3.94 X 10^3 


Ci 




492.161 


2.354 


2.23 


0.190 


0.412 


7.69 ± 1.45 X 10^^ 


cs 


2-1 


97.981 


0.294 


2.02 


0.306 


0.031 


3.38 ±0.45 X 10^2 


cs 


3-2 


146.969 


0.071 


1.47 


0.510 


0.013 


3.75 ±0.88 X 10" 


C34S 


3-2 


144.617 







0.207 


0.015 


< 1.13 X 10" 


CN 


1 3 3_Q 1 1 
' 2 ' 2 2 ' 2 


113.488 


0.032 


1.64 


0.660 


0.008 


2.29 ± 0.80 X 10^2 


CN 


1 3 5_Q I 3 
^' 2 ' 2 2 ' 2 


113.491 


0.042 


2.18 


0.660 


0.008 




SO 


Ifl-Oi 


30.002 


0.095 


1.80 


0.390 


0.018 


1.88 ± 0.45 X lO-"^^ 


so 


32-2i 


99.300 


0.133 


1.95 


0.302 


0.016 


4.17 ± 0.63 X 10^^ 


DCO+ 


1-0 


72.039 






0.416 


0.018 


< 1.17 X lO-*^-"^ 


jjl3cQ+b 


1-0 


86.754 








0.346 


0.019 


< 8.27 X 10"' 


HCO+ 


1-0 


89.189 


0.351 


2.74 


0.840 


0.035 


1.45 ±0.20 X 10^2 


HCO+ 


3-2 


267.558 






0.350 


0.013 


< 1.79 X 10"' 


DCN'^ 


1-0 


72.415 






0.414 


0.148 


< 1.65 X 10^2 


Ri^CN 


1-0 


86.340 






0.347 


0.014 


^ 1 03 X 10^^ 


HCN 


1-0 


88.632 


0.274 


1.54 


0.338 


0.032 


1.10 ±0.30 X 10^2 


HNC 


1-0 


90.664 


0.083 


1.74 


0.331 


0.018 


3.48 ± 1.13 X 10" 


020*= 


1-0 


72.108 






0.416 


0.021 


< 3.32 X 10^2 


C2H 


1-0 


87.317 


0.199 


2.89 


0.858 


0.019 


5.15 ±0.64 X 1013 


C2H 




87.329 


0.085 


1.86 


0.858 


0.019 




N2H+ 


1-0 


93.174 






0.322 


0.020 


< 9.67 X 10"' 


H2C0 


li,o-li,i 


4.830 


-0.256 


3.2 


0.303 


0.055 


5.61 ± 1.20 X 10^2 


H2C0 


2i,i-2i,2 


14.488 


-0.021 


1.3 


1.515 


0.011 




H2C0 


2i,2-li,l 


140.840 


0.085 


2.48 


0.532 


0.022 




H2C0 


2i,i-li,o 


150.498 


0.054 


1.83 


0.498 


0.008 




HDCO'" 


2i,i-li,o 


134.285 






0.223 


0.012 


< 6.44 X lO^" 


NH3 


1-1 


23.694 


0.040 


2.0 


0.500 


0.006 


5.77 ±0.87 X IO12 


NH3 


2-2 


23.723 


0.015 


2.6 


0.500 


0.005 




HC3N 


9-8 


81.881 






0.366 


0.021 


< 2.62 X 10" 


C3H2 


2l,2-lo,l 


85.338 


0.098 


1.92 


0.351 


0.039 




CH3OH 


2_i-l_iE 


96.739 


0.060 


0.93 


0.310 


0.012 




CH3OH 


2o-lo A+ 


96.741 


0.043 


1.49 


0.310 


0.012 




''Data taken with the ARO 12 m in 
''Data taken with the ARO 12 m in 


Jan. 2006. 
Feb. 2005. 
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Table 2 

Comparison of line temperatures and column densities for EC2 derived from Large 
Velocity Gradient (LVG) models and observed deconvolved integrated intensities (see 

section O WHERE TOTAL COLUMN DENSITIES ARE CALCULATED FROM INDIVIDUAL OBSERVED 

transitions). 



Molecule and 
Transition 




^mb 

(K) 


Tivg 
(K) 


Afmb(X) 


Mvg/ 

A^mb 




(X) 


nivg 


(X) 






n(H2) 


r 


Tex 

(K) 


CO 


1-0 




7.863 


7.850 


2.48 


X 


10i« 


1.0 


2.50 


X 


10" 


3.00 


X 


10- 


■3 


3, 


.20 


X 


10^ 


0.622 


20.40 


"CO 


1-0 




1.499 


1.500 


4.24 


X 


10" 


0.7 


2.76 


X 


10" 


4.00 


X 


10- 


-4 


3, 


.20 


X 


10^ 


0.056 


30.90 


CiSO 


1-0 




0.119 


0.124 


3.11 


X 


1014 


0.6 


1.96 


X 


loi-i 


3.10 


X 


10- 


-5 


3, 


.20 


X 


10^ 


0.004 


37.50 


CO 


2-1 




6.289 


7.720 


8.81 


X 


10" 


2.9 


2.55 


X 


10" 


3.00 


X 


10- 


3 


3, 


,20 


X 


10^ 


2.220 


13.60 


"CO 


2-1 




1.897 


1.810 


1.90 


X 


10" 


1.2 


2.31 


X 


10" 


4.00 


X 


10- 


-4 


3, 


.20 


X 


10^ 


0.388 


10.20 


CiSO 


2-1 




0.167 


0.157 


1.33 


X 


loi-i 


1.1 


1.42 


X 


1014 


3.10 


X 


10- 


■5 


3, 


.20 


X 


10^ 


0.033 


9.42 


CO 


3-2 




3.715 


3.760 


4.92 


X 


10" 


4.8 


2.37 


X 


10" 


3.00 


X 


10- 


■3 


3 


.20 


X 


10^ 


2.120 


10.50 


"CO 


3-2 




0.716 


0.699 


6.37 


X 


1014 


3.3 


2.09 


X 


10" 


4.00 


X 


10- 


-4 


3, 


.20 


X 


10^ 


0.271 


8.62 


c"o 


3-2 






0.056 


3.94 


X 


10" 


4.9 


1.91 


X 


1014 


3.10 


X 


10- 


■5 


3, 


.20 


X 


10^ 


0.020 


8.44 


CS 


2-1 




0.294 


0.297 


3.38 


X 


10" 


6.8 


2.31 


X 


10" 


3.70 


X 


10- 


■6 


1, 


.20 


X 


104 


1.370 


3.22 


CS 


3-2 




0.071 


0.071 


3.75 


X 


lOii 


45 


1.68 


X 


10" 


3.70 


X 


10- 


■6 


1, 


.20 


X 


104 


5.620 


3.00 


so 


lo-O 


1 


0.095 


0.096 


1.88 


X 


10" 


0.2 


3.11 


X 


10" 


5.60 


X 


10- 


■7 


1, 


.40 


X 


10^ 


-0.017 


-2.83 


so 


32-2 


1 


0.133 


0.129 


4.17 


X 


10" 


0.8 


3.37 


X 


10" 


5.60 


X 


10- 


-7 


1, 


.40 


X 


10^ 


0.045 


6.62 


HCO+ 


1-0 




0.351 


0.346 


1.45 


X 


10" 


1.9 


2.70 


X 


10" 


3.20 


X 


10- 


-7 


1, 


.20 


X 


104 


0.659 


3.57 


HON 


1-0 




0.274 


0.274 


1.10 


X 


10" 


0.9 


9.50 


X 


10" 


2.00 


X 


10- 


-7 


1, 


.40 


X 


10^ 


0.079 


4.50 


HNC 


1-0 




0.083 


0.080 


3.48 


X 


10" 


7.7 


2.68 


X 


10" 


5.00 


X 


10- 


-7 


1, 


.20 


X 


104 


0.656 


2.93 


N2H+ 


1-0 






0.028 


9.67 


X 


10" 


1.8 


1.70 


X 


10" 


2.75 


X 


10- 


-8 


1, 


.20 


X 


104 


0.142 


3.51 


H2CO 


li.o- 


-li.i 


-0.256 


-0.250 


5.61 


X 


10" 


2.5 


1.38 


X 


10" 


1.40 


X 


10- 


-6 


1, 


.20 


X 


104 


0.152 


0.95 


H2CO 


2i,i- 


-2l.2 


-0.021 


-0.021 








1.0 


5.61 


X 


10" 


1.40 


X 


10- 


6 


1, 


.20 


X 


104 


0.016 


1.43 


H2CO 


2i,2- 


-ll.l 


0.085 


0.096 








1.9 


1.07 


X 


10" 


1.40 


X 


10- 


■6 


1, 


.20 


X 


104 


0.623 


3.05 


H2CO 


2i,i- 


-li.o 


0.054 


0.049 








1.4 


7.90 


X 


10" 


1.40 


X 


10- 


■6 


1, 


.20 


X 


104 


0.501 


2.94 



Note. — LVG eolumn densities calculated with Afivg(X) = nivg(X) X 1 pc X AV. 
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Table 3 

Comparison of molecular abundances ratios relative t o HCQ+ and in EC2 (position A] 
L134N (center), TMC-1 (average) and L134N (range) (Idickens et AL.II2OO0I : IPratap et al" 
19971 ). Trans(lucent) cloud observations by (|TurnerI|2000| ). For EC2 T.. = 20 K, ^(Ha) 



1.2 X lO"' CM~^ AND A^(H2) = 7.4 X 10^^ CM~2 assumed in calculation of fractional abundances 

FROM individual OBSERVED TRANSITIONS (SEE SECTION ^ . 



A/rr%l or*! 1 1 








T,1 "^ztN 






' 1 V CI n Q 
±1 diio 




X/H2 




CO 


1- 


-0 


17078 


11000*^ 


7800^ 








3.3 X 


10- 


07 


"CO 


1- 


-0 


2919 


172^^ 


122^^ 


111-188 






5.7 X 


10- 


08 


c"o 


1- 


-0 


214 


22.00 


15.60 


14.2-24 






4.2 X 


10" 


09 


c"o 


1- 


-0 


<38.8 














<7.6 X 


10" 


10 


CO 


2- 


-1 


6070 


11000'^ 


7800*^ 








1.2 X 


10" 


07 


"CO 


2- 


-1 


1308 


172^ 


122^^ 








2.6 X 


10" 


08 


c"o 


2- 


-1 


92 


22.00 


15.60 








1.8 X 


10- 


09 


CO 


3- 


-2 


3389 


11000'^ 


7800*^ 








6.6 X 


10- 


08 


"CO 


3- 


-2 


439 


172^^ 


122^^ 








8.6 X 


10" 


09 


c"o 


3- 


-2 


<27.2 


22.00 


15.60 








<5.3 X 


10- 


10 


Ci 






53020 













1.0 X 


10- 


06 


cs 


2- 


-1 


2.33 


0.124 


0.320 


0.069-0.138 


1.1 X 10" 


08 


4.6 X 


10" 


11 


cs 


3- 


-2 


0.26 


0.124 


0.320 








5.1 X 


10" 


12 


c^-^s 


3- 


-2 


<0.08 












<1.5 X 


10" 


12 


CN'i 






1.58 


0.061 


0.070 


<0.045-0.069 






3.1 X 


10" 


11 


so 


1- 


-0 


13.0 


0.719 


0.130 


0.264-0.738 


3.2 X 10- 


08 


2.5 X 


10- 


10 


so 


3- 


-2 


2.87 


0.719 


0.130 


0.264-0.738 






5.6 X 


10- 


11 


DC0+ 


1- 


-0 


<0.08 












<1.6 X 


10" 


12 


H"C0+ 


1- 


-0 


<0.06 


C 


C 








<1.1 X 


10" 


12 


HCO+ 


1- 


-0 










2 X 10" 


09 


2.0 X 


10" 


11 


DCN 


1- 


-0 


<1.13 












<2.2 X 


10" 


11 


H"CN 


1- 


-0 


<0.07 


C 


c 








<1.4 X 


10" 


12 


HCN 


1- 


-0 


0.76 


0.925 


0.490 


0.555-0.968 


3.6 X 10" 


08 


1.5 X 


10" 


11 


HNC 


1- 


-0 


0.24 


3.251 


1.680 


1.324-3.963 


2.5 X 10" 


09 


4.7 X 


10" 


12 


Caff 


1- 


-0 


35.5 


0.288 


0.300 


0.171-0.333 


6.6 X 10" 


08 


7.0 X 


10" 


10 


N2H+ 


1- 


-0 


<0.07 


0.077 


0.013 


0.031-0.077 


- 1 X 10" 


09 


<1.3 X 


10- 


12 


H2CO 






3.87 








6.3 X 10" 


09 


7.6 X 


10" 


11 


NH3 






4.0 


7.635 


2.770 


4.284-9.127 


2.1 X 10" 


08 


7.8 X 


10- 


11 


HC3N 


9- 


-8 


<0.18 


0.054 


0.150 


0.030-0.073 


5 X 10" 


^0 


<3.5 X 


10- 


12 


CH3OH 


2- 


-1 




0.641 


0.099 


0.311-0.641 


1.8 X 10" 


08 









''A "CO/C^**0 ratio of 7.81 was assumed in L134N and TMC-1. 
''A ^^CO/C^^O ratio of 500 was assumed in L134N and TMC-1. 

i2q/13C ratio of 64 was assumed in L134N and TMC-1. 
'^From the sum of the two 113.49 GHz components assuming a relative intensity for this hne of 0.456. 
"From the 87.317 GHz component assuming a relative intensity for this line of 0.4167. 
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Table 4 

Range of dust to gas ratios and extinction for EC2 based on D and f3. 







1 


1.25 


1.5 


1.75 


2 


Q1200 




7.81e-5 


4.44e-5 


2.52C-5 


1.43e-5 


8.14e-6 






1.953 


1.110 


0.630 


0.358 


0.203 




= 14kpc) 


0.0015 


0.0027 


0.0048 


0.0084 


0.0148 


Mdust/Mgas {D 


20kpc) 


0.0011 


0.0019 


0.0033 


0.0059 


0.0104 


Aw (mag) (i?v 


= 3.1) 


4.6 


8.0 


14.1 


24.9 


43.8 


A\ (mag) (i?v 


= 2.0) 


2.9 


5.2 


9.1 


16.1 


28.3 



Note. — -Q^, an d k„ calculated from the 1983 Chicago Assumptions of 



Hildebrandl (|1983D 



Table 5 

Initial abundances relative to H. Most (C+, N, O, S, Si and Fe+) have been reduced by a 

FACTOR OF 5 FROM THOSE TYPICALLY USED TO MODEL THE LOCAL ISM. 



Species Abundance Species Abundance 



He 


1.4 X 10- 


1 


HD 


3.2 X 10 


-5 


C+ 


1.46 X 10" 


5 


N 


4.28 X 10 


-6 





3.51 X 10" 


5 


S 


2.0 X 10 


-8 


Si 


3.99 X 10" 


9 


Fe+ 


1.99 X 10 


-9 



Table 6 

ec2 chemical model agreement results ranked by fit as described in section [5l where 
n(H2) = 1.2 X 10** CM"^, T = 20K, in mag, CRI in terms of the standard ISM rate of 

1.3 X 10~^''S~-^, UV FIELD IN terms OF THE LOCAL ISM, AND MOST INITIAL ABUNDANCES (lA) 
REDUCED BY FACTOR INDICATED FROM TYPICAL LOCAL ISM VALUES. 



Model 


Fit 


Ay 


CRI X 


UV X 


lA X 


1 


0.01 


10 


1 


1 


0.2 


5 


0.11 


1 


20 


0.1 


0.2 


3 


0.16 


2 


10 


1 


0.2 


8 


0.40 


2 


20 


1 


0.1 


10 


0.41 


1 


20 


1 


1 


9 


0.42 


1 


20 


1 


0.4 


2 


0.43 


1 


20 


1 


0.2 


4 


0.43 


2 


10 


10 


0.2 


6 


0.43 


1 


20 


20 


0.2 


7 


0.44 


3 


20 


40 


0.2 
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Table 7 

Comparison of observations of EC2 with four different model predictions at 

STEADY-STATE. ABUNDANCES GIVEN RELATIVE TO A^(H2) = 7.4 X 10^^ CM"^, WITH n(H2) = 

1.2 X lO'^CM"^ AND T = 20K. Note that the abundances of the ^^C species have been 

MULTIPLIED BY 60 AND THAT OF THE C-^*0 SPECIES BY 500. MOST MODEL INITIAL ABUNDANCES 
RELATIVE TO H REDUCED BY A FACTOR OF 5 FROM TYPICAL LOCAL ISM VALUES. BOLD TYPE 
INDICATES AGREEMENT TO WITHIN A FACTOR OF 5. 



Species 


Observed 


iviociei i 


ivioaei z 


iVLociei 


iviodei 4 




m tii\^z 


Ay = 10 mag 


Ay = 1 mag 


Ay = 2 mag 


Ay = 2 mag 






L-xti — i X 


r~^T?T on \y 


UJXi — iU X 


v^xvl — iU X 






U V = i X 


U V = i X 


U V = i X 


U V = iU X 


"CO 


1.45e-06 


2.90e-05 


2.51e-07 


2.22e-06 


2.28e-07 


c 


1.04e-06 


1.41e-08 


1.37e-05 


2.47e-05 


1.79e-05 


c"o 


1.36e-06 


2.90e-05 


2.51e-07 


2.22e-06 


2.28e-07 


C2D 


<4.49e-ll 


1.14e-10 


1.48e-ll 


6.70e-ll 


1.18e-ll 


C2H 


6.96e-10 


1.65e-09 


1.71e-08 


2.79e-08 


1.04e-08 


CN 


3.09e-ll 


7.61e-09 


1.24e-09 


6.82e-09 


1.51e-09 


CO 


1.22e-07 


2.90e-05 


2.51e-07 


2.22e-06 


2.28e-07 


CS (2-1) 


4.57e-ll 


6.97e-09 


2.87e-10 


2.11e-09 


1.46e-10 


CS (3-2) 


5.07e-12 


6.97e-09 


2.87e-10 


2.11e-09 


1.46e-10 


DCN 


<2.23c-ll 


9.37e-10 


1.24e-14 


2.98e-13 


2.08e-14 


DCO+ 


<1.58c-12 


1.78e-09 


6.16e-15 


5.12e-13 


7.53e-15 


H"CN 


<8.35c-ll 


6.13e-09 


1.07e-ll 


1.28C-10 


1.36e-ll 


H13C0+ 


<6.71e-ll 


7.26e-09 


6.60e-12 


5.36e-ll 


5.16e-12 


H2CO 


7.58e-ll 


1.60e-08 


7.70e-ll 


1.19e-09 


9.38e-ll 


HC3N 


<3.54e-12 


1.55e-10 


7.26e-13 


1.19e-ll 


2.57e-13 


HCN 


1.49e-ll 


6.13e-09 


1.07e-ll 


1.28C-10 


1.36e-ll 


HCO+ 


1.96e-ll 


7.26e-09 


6.60e-12 


5.36e-ll 


5.16e-12 


HDCO 


<8.70e-13 


4.08e-09 


2.68e-13 


2.59e-ll 


4.44e-13 


HNC 


4.70e-12 


1.58e-08 


2.09e-ll 


4.05e-10 


2.48e-ll 


N2H+ 


<1.31c-12 


4.21e-10 


1.88e-14 


1.55e-12 


4.75e-14 


NH3 


7.80e-ll 


1.21e-07 


2.02e-13 


5.14e-ll 


3.04e-13 


SO 


1.34C-10 


7.16e-09 


6.43e-16 


1.47e-14 


6.16e-16 



Note. — For Tables [3 [8] and [9l Models only consider main isotopes (with the excep- 
tion of H and D), so for ^^CO, H^'^CN and H^^CO+, models are compared to observed 
abundances x 60; and for C^®0, models are compared to observed abundance x 500. 
For multiple transition detections of ^^CO, ^^CO and C^^O, a a weighted mean of each 
isotopomer's observed column density is used, where a — Tnih/rmsmh- 
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Table 8 

Comparison of observations of EC2 with three different model predictions at 

STEADY-STATE, BUT WITH VARYING UV FIELD. ABUNDANCES GIVEN RELATIVE TO iV(H2) = 

7.4 X lO^^CM"^, WITH n(H2) = 1.2 x 10'' CM~^ and T = 20K. Note that the abundances of the 

^^C species have been multiplied by 60 AND THAT OF THE C^^O SPECIES BY 500. MOST MODEL 

initial abundances relative to h reduced by a factor of 5 from typical local ism values. 

Bold type indicates ac;reemext to within a factor of 5. 



k_/ V/ J. wkj 




Model 5 


Model 6 


Model 7 




in EC2 


A\7' = 1 TTlrlE*' 


A\r = 1 TTlrlE*' 








cm = 20 X 


cm = 20 X 


cm = 20 X 






UV = 0.1 X 


UV = 20 X 


UV = 40 X 


i^CO 


1.45e-06 


2.79e-06 


6.66e-09 


3.39e-07 


c 


1.04C-06 


2.23C-05 


1.55e-06 


2.50C-05 


CiSQ 


1.36e-06 


2.79e-06 


6.66e-09 


3.39e-07 


C2D 


<4.49e-ll 


1.17e-10 


1.03e-12 


5.90e-12 


C2H 


6.96C-10 


6.97C-08 


1.34e-09 


3.90C-09 


CN 


3.09e-ll 


1.24e-08 


3.98e-ll 


8.46e-10 


CO 


1.22e-07 


2.79e-06 


6.66e-09 


3.39e-07 


CS (2 1) 


4.57e-ll 


7.17C-09 


3.30C-12 


5.92e-ll 


CS (3-2) 


5.07e-12 


7.17e-09 


3.30e-12 


5.92e-ll 


DCN 


<2.23e-ll 


4.51e-13 


2.27e-15 


1.13e-14 


DCO+ 


<1.58c-12 


6.42e-13 


8.48e-16 


1.86e-14 


H13CN 


<8.35e-ll 


2.44e-10 


2.13e-12 


7.56e-12 


H13C0+ 


<6.71e-ll 


8.63e-ll 


1.13e-12 


9.94e-12 


H2CO 


7.58C-11 


1.14C-09 


4.09C-12 


8.30e-ll 


HC3N 


<3.54e-12 


5.24e-ll 


l.OOe-15 


1.09e-13 


HCN 


1.49e-ll 


2.44e-10 


2.13e-12 


7.56e-12 


HCO+ 


1.96C-11 


8.63e-ll 


1.13C-12 


9.94e-12 


HDCO 


<8.70e-13 


1.98e-ll 


8.56e-15 


6.28e-13 


HNC 


4.70e-12 


5.57e-10 


2.78e-12 


3.68e-ll 


N2H+ 


<1.31c-12 


1.92e-12 


2.87e-17 


6.52e-14 


NH3 


7.80e-ll 


6.18e-ll 


3.83e-15 


1.35e-12 


SO 


1.34e-10 


2.96e-14 


1.56e-17 


9.25e-16 
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Table 9 

Comparison of observations of EC2 with three different model predictions at 

STEADY-STATE, BUT WITH VARYING INITIAL ABUNDANCES. ABUNDANCES GIVEN RELATIVE TO iV(H2) = 

7.4 X lO^^CM"^, WITH n(H2) = 1.2 x lO^CM"^ and T = 20K. Note that the abundances of the 

^^C SPECIES HAVE BEEN MULTIPLIED BY 60 AND THAT OF THE C-^*^0 SPECIES BY 500. MOST MODEL 
INITIAL ABUNDANCES RELATIVE TO H REDUCED BY A FACTOR (lA SHOWN BELOW) FROM TYPICAL 
LOCAL ISM VALUES. BOLD TYPE INDICATES AGREEMENT TO WITHIN A FACTOR OF 5. 



Species 


Observed 


Model 8 


Model 9 


Model 10 




in EC2 


Ay = 2 mag 


Ay = 1 mag 


= 1 mag 






CRI = 20 X 


CRI = 20 X 


CRI = 20 X 






UV = 1 X 


UV = 1 X 


UV = 1 X 






TA _ n 1 y 

— U.J. A 




TA — 1 V 


13C0 


1.45e-06 


3.48e-07 


6.36e-07 


2.21e-06 


C 


1.04e-06 


1.30e-05 


3.05e-05 


8.90e-05 




1.36e-06 


3.48e-07 


6.36e-07 


2.21e-06 


C2D 


<4.49e-ll 


3.58e-12 


1.22e-ll 


7.92e-12 


C2H 


6.96e-10 


1.92e-09 


1.99e-08 


1.98e-08 


CN 


3.09e-ll 


7.89e-10 


3.69e-09 


8.57C-09 


CO 


1.22e-07 


3.48e-07 


6.36e-07 


2.21e-06 


CS (2-1) 


4.57e-ll 


6.44e-ll 


9.73e-10 


3.46e-09 


CS (3-2) 


5.07C-12 


6.44C-11 


9.73C-10 


3.46e-09 


DCN 


<2.23e-ll 


8.97e-15 


2.10e-14 


3.19e-14 


DCO+ 


<1.58e-12 


3.12e-14 


4.19e-15 


2.51e-15 


Ri^CN 


<8.35c-ll 


6.31e-12 


3.04e-ll 


7.78e-ll 


H13C0+ 


<6.71e-ll 


1.06e-ll 


6.18e-12 


5.80e-12 


H2CO 


7.58e-ll 


7.36e-ll 


7.73e-ll 


6.19e-ll 


HC3N 


<3.54c-12 


1.85e-13 


1.19e-12 


1.24e-12 


HCN 


1.49e-ll 


6.31e-12 


3.04e-ll 


7.78e-ll 


HCO+ 


1.96e-ll 


1.06e-ll 


6.18e-12 


5.80e-12 


HDCO 


<8.70c-13 


9.07e-13 


1.68e-13 


7.19e-14 


HNC 


4.70e-12 


5.08e-ll 


4.23e-ll 


8.88e-ll 


N2H+ 


<1.31e-12 


4.78e-14 


6.39e-14 


1.65e-13 


NH3 


7.80e-ll 


6.46C-12 


2.80C-13 


3.08C-13 


SO 


1.34e-10 


1.35e-15 


1.91e-15 


6.54e-15 
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Table 10 
Summary of EC2 properties. 



Property 


Value 


Position A 


a = 48™ 38. 5^ 5 = 58° 28' 28.1" (J2000) 




/ = 137.76, 6= -0.98 


Radial velocity 


Viad = —103.70 kms~^ 


Galactocentric distance 


i? = 22 to 28kpc 


Distance from Galatic plane 


250 to 350 pc 


Size 


6»EC2 = 30 to 40 pc 


Temperature 


T = 20K 


Density 


n(H2) - 10'* cm-^' 


Cosmic ray ionisation rate 


cm = 10-20 X 1.3 X 10-1^ s-i 


UV photon field 


UV = 10-20 X local ISM values 


Initial abundance 


lA = 20 per cent of local ISM values 


Extinction 


Ay < 4 mag 


Dust to gas ratio 


Mdust/Mgas = 0.001-0.015 


Mass 


Mec2 « 10^ Mq 
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-100 - 



-200 - 



-300 - 



-400 





Fig. 1.— JCMT 15m May- July 2005 combined 
contour map ol observed CO 2-1 intensities, show- 
ing observed positions A through K. Positions A, 
J and K correspond to pos itions copkl , copk 2 
and copkS respectively from Digel et al] (1994). 
Axis offsets in arcsec from position E (a20oo = 
02:48:38.5, S = 58:28:58.3, Wrad -103.70 
kms-i). EC2MAP1 (top) CO 2-1: peak = 
6.906 K, rms = 0.321 K. EC2MAP2 (bottom) CO 
2-1: peak = 7.124 K, rms = 0.22 K. 



EDGE2A (02 44 52.6/+58 16 00) 
23.723 GHz n NH3 (J,K)=(2,2) 




23.694 GHz 



NH3 (J,K)=(1,1) 



EDGE2A (02 44 52.6/+58 1 6 00) 

14.488 GHz 5 X H2C0 2(1 1)-2(12) 





Fig. 2.— EC2 spectra observed at the MPIfR Ef- 
felsberg 100 m position A: (top) NH3 (J,K) = (1,1) 
and (2,2) (HF = group of satellite hyperfine com- 
ponents), (middle) H2CO 2i,i-2i,2 and li,o^li,i- 
Note that the scale for the 2i_i-2i_2 transition has 
been multiplied by 5. (bottom) SO Iq^Oi. 
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HC0+PKEDGE2 


110.201353 GHz 




0.15 


12M-FB12 


Jin 13CO(1-0) 










t; (k) 












HC0+PKEDGE2 
12M-FB12 




-120 


Resl freqijency (UHz) 
8.919 10 e.91S5 10 


- EDGE2P 
12M-FB22 


89.188518 GHz 
pn HC0+(1-0) 





HC0+PKEDGE2 


1 13.490982 GHz 




12U-FB22 








f\ f" 


0.05 

















HCG+PKEDGE2 
12U-FB12 



99.299879 GHz 
SG(3 2-2 1) 





1.50505 10 1.505 1i 



Rejt Frsouency (MHz) 



HC0+PKEDGE2 
12M-FB22 



150.49839 GHz 
H2C0(2 1 1-1 1 0) 



HC0+PKEDGE2 
12U-FB12 



85.3382 GHz 
C3H2(2-1) 



Rest Frequency (MH;) 




Fig. 3.— EC2 spectra observed at the ARO 12 m Feb-Junc 2002 position A (a2ooo = 02:48:38.5, 
5 = 58:28:28.1, Vrad = — 103.70kms^^). The ordinate is T|j.(K), uncorrected for beam-fiUing. 
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Rb5| FrequEncy (UH;) Rest Frequency [MH^J 

4921B 10 4.9216 10 4.9214 IQ 2.6757 10 2.6756 10 2 6755 10 2.6754 IQ' 




Velocity (kmA) Velocity (Km/s) 



ig. 4.— EC2 spectra observed at the JCMT 15 m June 2004 position A (a20oo = 02:48:38.5, 
= 58:28:28.1, Wiad = -103.70 kms"!). The ordinate is T;^(K), uncorrected for beam-fiUing. 
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EC2MAP1 C0(2-l) 230.538 GHz JCMT-HAWAII 
Scans 514109-515107 

htensity: peak Tr = 6.906 K, rms = 0.321 K (contours 10%) 




EC2MAP1 00(3-2) 345.796 GHz JGMT-HAWAII 
Scans 523097-529050 

rtensity: peak = 4.457 K, rms = 0.75 K (contours 10%) 



0) 




-100 



50 





A RA 



-50 



EC2MAP1 1.5CO(7-1) 220.399 GHz JCMT-HAWAI 
Scans 515' ' 1-517066 

htensity: peak = .1033 K, rns = 0.765 K (contours 10%) 



0) - 




EC7MAP' 13CO(3-2) 330.588 GHz JCMT-HAWAII 
Scans 530058-629261 

htensity: peak = '.S79 K, rms = 0.12 K (contours 10%) 



0) 



-100 




Fig. 5. — JCMT 15 m May- July 2005 maps of observed CO intensities centered on position E (a2ooo = 
02:48:38.5, S = 58:28:58.3, u^ad = -103.70 kms"!) with position A (azooo = 02:48:38.5, S = 58:28:28.1) 30" 
below map center. Axis offsets in arcscc, contour levels 10% of peak, dashed contours 50% of peak (FWHP) 
and white circle is half-power beam width (HPBW). (a) CO 2-1: peak = 6.906 K, rms = 0.321 K. (b) 
CO 3-2: peak = 4.452K, rms = 0.25K. (c) ^^CO 2-1: peak = 3.033K, rms = 0.265K. (d) ^^CO 3-2: 
peak = 1.579 K, rms = 0.12K. 
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200 200 



Fig. 6.— IRAM 30 m Oct-Nov 2005 1.2 mm dust 
map centered on position EC2DUST1A (a20oo = 
02:48:41.0, 5 = 58:29:27.9) with CO 2-1 inten- 
sity contours from Fig. [5^ superimposed. The 
rms noise for the map is 6.35 mJy /beam with a 
detected peak of 20.3 mJy /beam. Axis offsets in 
arcsec. 




100 50 -50 -100 



Fig. 7. — CO 2-1 intensity contours from 
EC2MAP1 (Fig. [Hi) overlaid on 1.2 mm dust map 
(Fig. [6] EC2DUST1A), showing the correspon- 
dence of peak CO 2-1 and 1.2 mm dust emission. 
Axis offsets in arcsec. 



27 




Fig. 8. — Agreement factor for models with varying: (a) (mag) and cosmic ray ionisation (CRI); (b) 
Ay and CRI with UV field increased by a factor of 10 compared to local ISM values; (c) Ay and CRI with 
UV field reduced by a factor of 10; (d) Ay and an increasing UV field; (e) Ay and most initial abundances 
reducing relative to typical local ISM values; (f) Ay and n(H2). Unless indicated above, parameters were 
fixed at: T = 20K; ^(Ha) = 1.2 x lO^^cm^^. ^j^j ^.3,^ = 20x the standard ISM rate of 1.3 x lO^^^s-i; UV 
field set to local ISM values; most initial abundances reduced by a factor of 5 from typical local ISM values 
(see Table [5]). Numbered disks indicate specific models relative to the parameter surface. 
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- HCN 
NHa 





CeH 
HNC 



SO 
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Fig. 9. — Fractional abundances varying over time (see Tables [71 [8] and [Ql) for: 

(a) Model 2, where Ay = 1 mag, CRI = 20x , UV = 1 x , lA = 0.2x, n(H2) = 1.2 x 10'' cm^^^ T = 20K; 

(b) Model 4, where Ay = 2 mag, CRI = lOx, UV = lOx, lA = 0.2x, ^(Ha) = 1.2 x 10^ cm'^, T = 20K; 

(c) Model 6, where Ay = 1 mag, CRI = 20x, UV = 20x, lA = 0.2x, ^(Hs) = 1.2 x 10'' cm'^, T = 20K; 

(d) Model 7, where Ay = 3 mag, CRI = 20x, UV = 40x, lA = 0.2x, n(H2) = 1.2 x lO^cm^^^ T = 20K; 

(e) Model 8, where Ay = 2 mag, CRI = 20x, UV = Ix, lA = O.lx, n(H2) = 1.2 x lO^cm-^, T = 20K. 

(f) Agreement factor varying over time for Models 1-4, where Ay (mag) and cosmic ray ionisation (CRI) 
rates are varied (see Figs. [S^,b and Table [7]). 
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